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1. Introduction 

The use of "partially quenched" QCD (PQQCD), in which valence and sea quark masses 
are varied independently, has become common practice in large-scale computations of hadronic 
quantities from QCD using lattice methods. This usefulness arises from the following observations: 

1 . Sea quark masses are, by definition, the quark masses that go into the generation of dynam- 
ical gauge field configurations on which hadronic quantities are evaluated. Valence quark 
masses are, by definition, the quark masses that go into the operators that are evaluated on 
these gauge configurations. Such operators are contractions of quark propagators, which de- 
pend on the gauge field configurations. For many applications, it is numerically less expen- 
sive to generate quark propagators than gauge field configurations. Given limited resources 
it can thus be advantageous to generate data for many values of the valence quark masses, 
using configurations depending on only a limited number of sea quark mass values. 

2. PQQCD with a given number of sea quarks contains full QCD with the same set of sea 
quarks It follows that the low-energy constants (LECs) of the low-energy effective the- 
ory are those of the real world, because by construction, these LECs are independent of the 
quark masses [g]. For the case of unquenched QCD, this low-energy theory is chiral pertur- 
bation theory (ChPT), which provides a systematic expansion in (light) quark masses, with 
coefficients characterized by these LECs. 

3. On the lattice, PQQCD can be generalized to "mixed-action" QCD, in which not only va- 
lence and sea quark masses are independently chosen, but also the discretization of the Dirac 
operator is chosen independently in the sea and valence sectors of the theory. This is a gener- 
alization in the sense that the continuum limit of such a lattice theory is described by PQQCD 

In order to apply these ideas, one needs the correct effective theory for PQQCD, i.e., one 
needs the partially quenched version of ChPT (PQChPT) [Q]. While it is relatively straightforward 
to extend ChPT to the partially quenched case, it is less clear than in the unquenched case whether 
indeed PQChPT is the correct effective theory for PQQCD. This is because PQQCD violates a 
number of basic properties of a healthy quantum field theory on which arguments that ChPT is the 
correct effective theory for unquenched QCD rely. This follows from the definition of PQQCD as 
a euclidean path integral that includes an integral over ghost quarks, which have the same quantum 
numbers as the valence quarks, but the opposite (and thus "wrong") statistics. 

For many field theories, the path integral for a euclidean theory can be expressed in terms of 
a transfer matrix. If the field theory is healthy, this transfer matrix is hermitian and bounded, and 
contact with the hamiltonian formulation of the theory can be made. 2 But if the theory contains 
fields with the wrong statistics, it is not clear what properties of such a construction survive. 

In Ref. [g], the validity of ChPT as a systematic low-energy effective theory for the Goldstone 
sector of QCD was conjectured to follow from the basic properties of a healthy quantum field 

'They do depend on the number of sea quark flavors. 

2 Of course, quantum field theories exist which are believed to be healthy, but for which no transfer matrix can easily 
be constructed. 
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theory: analyticity, unitarity, cluster decomposition, and symmetry. Stated differently, the S-matrix 
calculated with the most general local lagrangian consistent with a certain symmetry group was 
conjectured to be the most general possible 5-matrix consistent with these basic properties. This 
was then used as a starting point for the systematic development of ChPT as an expansion of S- 
matrix elements in terms of the pion momenta, following a well-defined power-counting scheme. 
The reliance of this argument on unitarity, though, appears to be a fundamental difficulty in trying 
to extend it to the partially quenched case, which is certainly not unitary. 

An alternative justification for ChPT as the low-energy effective theory for QCD was presented 
in Ref. [Q]. The argument detailed there was based on locality and clustering of the underlying the- 
ory of unquenched QCD, as well as its symmetries. Locality and clustering guarantee the existence 
of vertices in the effective theory that are independent of the correlation functions in which they 
appear, and, consequently, the existence of a loop expansion. This approach seems more fruitful for 
PQQCD than that of Ref. [g]. By construction, PQQCD is local. It is less clear that it also satisfies 
the cluster property, but numerical evidence suggests that that is indeed the case, because euclidean 
correlation functions are observed to decay exponentially in distance (up to possible nonstandard 
powers, which come from the characteristic double poles |[l], [7|] in the theory). 

The first requirement for the chiral theory of PQQCD is, of course, the dynamical breaking 
of chiral symmetry in the theory with massless quarks. The Goldstone bosons associated with 
this breaking provide the low-energy degrees of freedom for the low-energy effective theory. It 
was argued in Ref. [|7|] that chiral symmetry must also be broken in PQQCD. The argument is 
that this must happen in the sea sector of that theory, because it is identical to unquenched QCD. 
This follows, since, by construction, the valence (and ghost) quarks are not part of the dynamics. 
Goldstone bosons corresponding to pions made only out of sea quarks thus have to be present 
in the partially quenched theory. Furthermore, the vectorlike partially quenched symmetries that 
relate sea and valence quarks then imply that there also have to be Goldstone excitations in the 
valence and ghost sectors, thus providing the necessary degrees of freedom for the construction of 
PQChPT. 3 

The key outstanding issue for the justification of PQChPT as the correct chiral theory for 
PQQCD therefore seems to be the question of whether PQQCD obeys the clustering property. Here 
we begin an investigation into the clustering property of PQQCD by considering the construction 
of a transfer matrix for PQQCD. We start with the ghost sector, which is the ultimate source for the 
nonstandard features of the partially quenched theory. 

2. Staggered ghosts 

PQQCD contains three types of quarks: sea quarks, valence quarks, and ghost quarks, which 
have the same masses as the valence quarks, but opposite statistics. Because of this quark content, 
the complete fermion determinant of the partially quenched theory is just that coming from the sea 
sector, while the valence and ghost determinants cancel each other. 

Since ghost quarks violate spin-statistics, we expect that obstructions to the existence of a 
well-behaved transfer matrix would originate from this sector of the theory. We therefore consider 

3 These vectorlike symmetries are not spontaneously broken fi|. 
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a theory of bosonic staggered quarks, in an arbitrary gauge field background. 4 If we choose 

S = J £xHx)D(x,y)x(y) + Z m xHx)x(x) (2.1) 

x,y x 

with 

r llx (x) = (-ir + - +x ^ , (2.2) 

then 

Z = jDx^Dxexpi-S) (2.3) 

is convergent for m > 0, because the staggered Dirac operator D is anti-hermitian. The field % is 
bosonic, i.e., it is a c-number valued staggered quark field. Of course, S exhibits the same species 
doubling as the normal staggered quark action, and therefore we employ a two time-slice method 



for constructing the transfer matrix representation of Z, guided by the construction in Ref. [10]. 

We split X into its real and imaginary parts by defining %(x) = Tj4(x)0i (x) + i<j>z(x), and then 
identify 

t = 2k : fafct) = 3>i,ft(x) , fafct) = -®2,k(x) , (2-4) 
t = 2k+ 1 : <j>i(x,t) = U 2: k(x) , fa(x,t) = Ui,k(x) . 

Then, if the extension of the (periodic) lattice in the time direction is T (which we take to be even), 
the partition function can be written as the trace of the T/2-th power of a transfer matrix, 

/T/2 \ 

Z(4r)=Trm?tWl > (2-5) 

in which 

t k (ty) = g-*i^-[^(2(fe+i)]*2 e -n 2 ^r + [^(2fe+i)]n 1 _ (2 6) 

Here J#± are defined through 



WiJ^lWit)}^ = ES ±I>K*) Re Uj(x,t)y 2 {x + j) + (1++ 2) J (2.7) 

- £ (x) (Vj (x) Im ^ (x, O^i (x + j) - ( 1 -> 2)' 



+m(^ 1 (x) 2 + ^ 2 (x) 2 )| , 
= *b(*)f?4(x) , 

and the hermitian operators 4>i 2 and TJ12 satisfy the commutation rules 

[& a (x),n b (y)] = iS(x-y)8 ab . (2.8) 



4 We use staggered quarks here because the nonperturbative definition of PQQCD with staggered quarks is straight- 
forward. For other types of quarks, see Ref. H and refs. therein, in particular Ref. Wm. 
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From f, we can define a hamiltonian in the limit in which the temporal lattice spacing goes to zero; 
with 

H[W(t)] = lim -logf*(^)/(2a t ) , (2.9) 

h, = im£(n2(x)+n2(x)+<i>2(x)+4>2(x)) , 

x 

H 2 [^(t)] = \y\ V'jW (n 2 (x)Ref/X^0ni^+7)-<J>i(^)Ref/,(^04>2^+J)) 
-r\j{x) (-fli(x) Im U j (x,t)tl l (x + j) + 4>i(x) Im Uj{x,t)&x (x + jfj 
+(lo2)| . 

Both H\ and H2 are hermitian, and do not commute. Therefore, H is not hermitian, and not normal. 
It follows that the transfer matrix is not hermitian and not positive definite. 

Despite these (nonsurprising) conclusions, it is possible to make progress. Our transfer matrix 
factorizes as 

T k (W) = T l {W)T 2 (W) , (2.10) 

in which T\ is the exponential operator with M 1 - and Ti is that with Both T\ and T2 are normal 
and bounded. Since the hamiltonians JffL and M' + have a positive real part (the parts proportional 
to m): 

||7l 2 ||<1, (2.11) 

which implies 

II t \\<\\ fi (HI f 2 II < 1 , (2.12) 

which establishes that all eigenvalues of f have an absolute value smaller than one. 5 It follows 
that correlation functions in this theory decay exponentially with distance if the eigenvalue Ao with 
maximal |Ao| is unique. Of course, the construction needs to be extended to the complete partially 
quenched theory, but we do not expect any difficulties with the other building blocks of a complete 
transfer matrix. 



3. Trivial gauge field background 

For a trivial gauge field background, it is straightforward to diagonalize H. In terms of canon- 
ical creation and annihilation operators, introduced through 

*.(.«) = / 75 (<■!(*)+<■!(-*)) e* 1 , (3.1) 

n 'W = /(0 ^ («.<*) - «!<-*)) «** • 

5 This product inequality is satisfied if we use the euclidean norm, which, for a matrix A, is defined as the square-root 
of the largest eigenvalue of A* A. 
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n 2 (^) = / ^3 ^ (a2(-t) + «4(*)) , (3.2) 



we find that 



^ = /| m ( a i^ ai ^ + 4(^) a 2(p)) + *X sin (/?,■) ^a 1 (p)a- / a 2 (p)-4(p)a ; a]'(p)^ j , (3.3) 

with —7i/2<pj<7t/2. The a ; are anti-hermitian 8x8 Dirac matrices because there are 2 3 spatial 
doublers. (The two doublers in in the time direction appear explicitly.) Now £ysin(py)a ; can be 
diagonalized, with eigenvalues 



±is(p) = ±i /£ sin 2 (pj), (3.4) 

making H a sum of terms of the form 

h(p) = m(a\(p)ai(p)+al(p)a 2 {p)) ±s(p) (ai(p)a 2 (p) -a\{p)a\{p) ) . (3.5) 



For each p, this can be diagonalized with a generalized Bogoliubov transformation [[TTJ] : 

b\ = cos0ai — sinGal > (3-6) 



2 ' 
t 



&2 = cos — sin a J 
&i = cos a] + sin a 2 , 
t>2 = cosQa^ + sinflai , 
where we note that hi is not the hermitian conjugate of b{. With 6 = ltan _1 (j/m) this yields 



h = E(b 1 bi+b 2 b 2 ) , E = Vm 2 + s 2 . (3.7) 

The operators bj and bj are annihilation and creation operators, and indeed, H has complete sets of 
left and right eigenstates. The form of the eigenvalues E are not a surprise: the ghost determinant 
should cancel the valence determinant, so we expect the eigenvalues for the ghost hamiltonian to 
match those of the valence hamiltonian. 

We can also calculate correlation functions in the theory with a trivial background. The two- 
point functions in a theory with hamiltonian (3.3) are 



(aj(t)aj{0)) = Sij^re- 

( ai (t)aj(0)) = -<«J(0«J(0)> = S i+j ,3^e- Et . 

While these correlation functions exhibit the expected exponential decay, they also clearly show a 
violation of unitarity. In a healthy field theory, by inserting a complete set of states between the 
operator at time t and the one at time 0, one would conclude that (aJ(t)aj(O)) has to be nonnegative. 
Here this is not the case, because the vacuum bra (0| defined by (0\bj = is not the conjugate 
of the vacuum ket |0) defined by bi\0) = since bj is not the hermitian conjugate of bj. The 
results (3.8) also follow directly from a path integral for this hamiltonian, without first performing 
the Bogoliubov transformation. 
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4. Discussion 

In this talk, we took the first, and most important, step toward the construction of a transfer 
matrix for PQQCD. We considered PQQCD with staggered fermions, because it is straightforward 
to define the ghost sector in that case. The ghost quarks are then "staggered ghosts," i.e., staggered 
quarks with bosonic statistics, and we constructed the transfer matrix for staggered ghosts in an 
arbitrary gauge field background. The rest of the construction remains to be done, but we anticipate 
that adding the sea and valence parts, as well as the gauge part, to the transfer matrix will be 
straightforward. 

We then proved that the ghost transfer matrix in an arbitrary gauge field background is bounded 
in the sense that the absolute value of its largest eigenvalue is bounded by one. At nonzero quark 
mass, correlation functions of ghost quarks are thus expected to decay exponentially, even if they 
do not satisfy the positivity properties of a healthy quantum field theory. We would expect this to 
carry over to PQQCD as long as that theory has a nonvanishing mass gap. In the case of PQQCD, 
the meaning of this would be that there is a unique eigenvalue Ao of the complete transfer matrix of 
the partially quenched theory with |Ao| maximal, while all other eigenvalues A satisfy |A| < |Aq|. 

We then checked this in the free ghost theory, in the limit of vanishing temporal lattice spacing. 
Indeed, we find that there is a nonzero mass gap, set by the quark mass m, and all two-point 
functions decay exponentially. This supports the conjecture that the argument of Ref. [Q] for the 
validity of ChPT can be carried over to the partially quenched case. Already in the free theory we 
find violations of unitarity in two-point functions, so it is clear that also PQQCD, and therefore 
PQChPT, will suffer from the same disease. However, PQChPT would still be the correct effective 
theory for PQQCD at low energy, because the underlying theory is local, and still satisfies cluster 
decomposition, despite its nonhermitian transfer matrix. 

We thank Michael Ogilvie for helpful discussions. 
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